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Abstract. We have undertaken a differential photometric 
I band search for variability in three Pleiades brown dwarfs 
and three very low mass field L dwarfs. Analysis of the re- 
sultant time series for the Pleiades brown dwarfs (Teide 1, 
Calar 3, Roque 11) gives no evidence for variability (above 
99% confidence) greater than 0.05 mags in any of these ob- 
jects on timescales between 25 minutes and 27 hours. Time 
series of two of the L dwarfs, 2MASSW J0913032+184150 
and 2MASSW Jl 146345-1-223053, are also consistent with 
no intrinsic variability (no significant variability above 
0.07 mags for 30 mins < r < 126 hours, and above 0.025 
mags for 30 mins < r < 75 hours, for the two objects 
respectively). However, the L dwarf 2MASSW J1145572 
-1-231730 shows evidence for variability of amplitude 0.04 
magnitudes with a tentative period of about 7 hours. It 
is therefore interesting that Kirkpatrick et al. (199£) re- 
port an Ha emission line of 4.2A equivalent width for this 
star: If our detection is confirmed as a rotation period, it 
would support the link between rotation and Ha emission 
in very low mass stars. 

Key vifords: methods: observational - stars: activity 
- stars: low-mass, brown dwarfs - stars: individual: 
2MASSW J1145572-h231730 - stars: rotation 



1. Introduction 

Brown dwarfs are of great interest to astronomers. Ac- 
cording to one definition they can be considered to bridge 
the gap between the lowest mass stable-hydrogen-burning 
stars and the Jupiter-like planets which do not get hot 
enough to burn deuterium, the most easily "combusted" 
element. Based on this definition, current models place 
their mass in the range w 0.015-0.0751/© (DAntona & 
Mazzitelli 1994; Burrows et al. 1997). However, for sin- 



gle objects mass is not a directly observable quantity, so 
brown dwarf status can only be conferred via secondary in- 
dicators. Models predict that brown dwarfs spend most of 



their lives cooling, and all but the youngest and most mas- 
sive brown dwarfs have effective temperatures below about 
2800 K. Of particular interest, therefore, are the recently 
discovered compact objects with atmospheres character- 
istic of very low temperatures (Kirkpatrick et al. 1999). 



Models of these objects suggest 1500 K < T^s < 2000 K, 
and the newly introduced class L dwarfs has been sug- 
gested to emphasise the status of these objects as a low 
temperature continuation of the M dwarfs. We therefore 
expect many of them to be brown dwarfs. 

The intrinsic properties of brown dwarfs and L dwarfs 
are of astrophysical interest, and the observations and 
models to date raise a number of important questions. For 
example, what is the nature of their atmospheres? At tem- 
peratures below about 2500 K, models predict that dust 
grains start to condense and appear in the atmosphere 



(e.g. AUard et al. 1997). How does the angular momen- 



tum of these objects evolve? Knowledge of their rotation 
speeds as a function of age and mass will help address 
the problem of the formation and evolution of these ob- 
jects. Do these objects have magnetic activity and chro- 
mospheres? In higher mass objects, dynamo models give 
a relationship between rotation speed and magnetic activ- 
ity (as observed, for example, via Ha emission). The mea- 
surement of rotation periods and variability amplitudes is 
important for addressing these questions. 

To date there has been very little work on measuring 
variability in brown dwarfs and L dwarfs. Tinney & Tol- 



ley (1999) provide some evidence for photometric variabil- 
ity in the M9 brown dwarf LP 944-20, but do not derive 
a period for what they suspect to be a rotational modu- 



(1997) 



for the L5 field brown 
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lation. Based on spectroscopy, Tinney & Reid (199S) re- 
port wsini of about 28 ± 2kms^^ for this object. Martin 
et al. (1998) give a vsini of 13 ± lkms~^ for the M6 
Pleiades brown dwarf Teide 2, and Martin et al 
derive a wsini of 20 ± lOkms 
dwarf DENIS-P J1228. 2-1547. 

In this paper we present the results of a search for vari- 
ability in six very low mass (VLM) stars (Table ll|) . Three 
of these are Pleiades brown dwarfs. Two (Teide 1 and 
Calar 3) are confirmed brown dwarfs, based on their radial 
velocities and lithium detections. The third, Roque 11, is 
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Table 1. Properties of brown dwarf and L dwarf targets. All data for the 2MASS objects 
are from Kirkpatrick et al. (1999). For these objects / is their Igpec value (magnitude 
evaluated from a flux-calibrated spectrum). 



lAU name 



SpT Ha EW 



LiiA6708 EW 

A 



2M0913 
2M1145 
2M1146 
Roque 11" 
Teide r 
Calar 3-'^ 



2MASSW J0913032+184150 
2MASSW J1145572-H231730 
2MASSW J1146345-H223053 
RPL J034712-f-2428.5'' 
TPL J034718-(- 2422.5'' 



< 1.0 

< 0.4 
5.1±0.2 

l.OiO.2-'^ 
18.73» Mg'' 6.5^-10.2'* ±1.0 1.8±0.4-'^ 



.07 


L3 


< 0.8 


.62 


LI. 5 


4.2 ±0.2 


.62 


L3 


<0.3 


.75'' 


M8" 


5.8 ±1.0" 


.80'' 


MS'' 


3.5-8.6'= 




"Zapatero Osorio et al. 

''Zapatero Osor io et al. 

''Rebolo et al. (|l99J) 

''Martin et al. ( [1996[ ) 

*= several values are given in the literature: 6.1 ± 1.0 (Rebolo et al. 1995); 



4.5 ±1.0 (Reb olo e t al. |1996| ); 3.5 ±2.0, 3.7 ±2.0, 8.6 ±2.0 (Zapatero Osorio et al. |l997a| ) 
^Rebolo et al. ( |l99e| ) 
''Zapatero Osorio et al. ( |l997b| ) 



a probable brown dwarf and cluster member. The three 
other objects (2M0913, 2M1145 and 2M1146) are L dwarfs 
discovered by the 2MASS survey. One of these, 2M1146, 
is a brown dwarf based on the detection of lithium in its 
spectrum. The other two have no or very little lithium, 
so are probably either VLM stars (just above the hydro- 
gen burning limit), or old high mass (> O.OGSM©) brown 
dwarfs. 

The most likely cause of variability would be rotational 
modulation of the emitted flux. In the case of the three L 
dwarfs, this could be due to inhomogenous dust or other 
clouds rotating across the unresolved stellar disk. In the- 
ory, a modulation could also be observed on the timescale 
of formation and dissipation of the clouds, although this 
may be longer than the few days duration of our search. 
For all objects, modulation could be caused by surface star 
spots induced by a magnetic field. 

2. Observations and Processing 

The project was carried out with the 2.2m telescope at 
Calar Alto, Spain, using CAFOS and the IKxlK TEK13c 
CCD detector. The observations were obtained over the 
six nights 08/09 January to 13/14 January 1999 (MJD 
2451187.4266 to 2451192.7560). However, due to poor 
weather (including frequent fog and poor seeing) and vari- 
ous technical problems, less than half of the available time 
could be spent obtaining data. The observing procedure 
was to cycle around the science fields in turn, taking five 
minute exposures in the I band in each field (Teide 1 and 
Roque 11 could be observed in the same field). A num- 
ber of other fields were included for calibration purposes. 
The spacing between observations of a given field was non- 
uniform, although was typically 15-30 minutes on a given 
night. 



The CCD frames were processed using IRAF. To help 
remove fringing in the CCD, a night sky flat was con- 
structed from a clipped combination of 31 frames of 15 
different fields. 

For each science field, reference stars were selected 
which satisfied the following conditions: 

— have a near-Gaussian, near-circular PSF 

— are bright, but never saturated 

— are isolated from other objects 

— are as near as possible to the science object. 

The reference stars were not know a priori to lack vari- 
ability. About ten reference stars were typically found in 
each field. Fluxes at a range of aperture sizes were then 
evaluated for each reference star and science star using 
the IRAF aperture photometry task phot. Sky subtraction 
assumed a single value for the sky across the extraction 
aperture. 

During reduction, 2M1146 was suspected to be a bi- 
nary on account of the shape of its PSF. This was con- 
firmed by Kirkpatrick et al. (1999) who report a back- 
ground star of much earlier spectral type separated by 
w 1" from 2M1146. 



3. Time Series Analysis 

Intrinsic variations in a science star were searched for by 
looking for changes in the differential magnitude of the 
science star, that is, changes in the magnitude of the sci- 
ence star relative to a number of reference stars. This 
makes the observations relatively insensitive to changes in 
atmospheric transparency during or between exposures. 
For this to be successful, we must be confldent that the 
reference stars are insignificantly variable themselves. We 
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therefore adopted the following procedure to test for vari- 
ability in the reference stars. A reference star, s, is chosen, 
and the average flux of all other reference stars formed. 
The relative magnitude of s is then evaluated with respect 
to this average. This is done at every epoch, producing a 
time series rus i, nis.2, • • • , "m-s t, ■ ■ -, from which the mean 



is subtracted so that J2t 



l'S,t 



0. At each epoch the 



total uncertainty in the relative magnitude, Cg,*, is also 
estimated (see below). We then use a x^ test to decide 
whether the scatter of these relative magnitudes is con- 
sistent with their measurement uncertainties, using the 
statistic 

2 



Xs 



t 



es,t 



(1) 



The null hypothesis for this test is that there is no vari- 
ability in the time series. If P(Xs) is the probability that 
the null hypothesis is true, then the larger xi the smaller 
P(x^). A star will fail the test (show evidence for vari- 
ability) if F{xi) < 0.01, a 99% confidence interval. (This 
would correspond to 2.5a significance if the distribution 
were Gaussian.) xi is measured for every reference star, 
and if any fails the test, the star with the highest x^ is 
rejected from the set of reference stars. This whole process 
(starting from forming the time series for each star) is it- 
erated until all remaining reference stars pass the x^ test. 
These then define the reference system from which to de- 
termine the relative magnitude time series of the science 
star. 

The point of using more than one reference star to 
determine the time series for the science star is that it 
increases the photometric precision of the relative magni- 
tudes. Furthermore, any small residual variability in any 
one reference star will have a smaller effect on the average. 

The x^ test is a good initial test for variability, as it 
is sensitive to any type of variability and not just peri- 
odic variations. However, for the test to be successful, the 
errors, es^t, must be accurately determined. Using photon 
statistics we can evaluate the "formal" errors in any mag- 
nitude determination. These are (1) photon arrival noise 
from the object, (2) photon arrival noise from the sky, and 
(3) imperfect determination of the subtracted sky level. 
There are, however, additional "informal" errors which 
need to be accounted for if we are not to overestimate x^ ■ 
The most important of these results from imperfect flat 
fielding. We have estimated (from the construction of dif- 
ferent flats) the accuracy of the flat to be 0.5%, and have 
combined this in quadrature with the formal error sources 
to give Cs^t- Other informal error sources include (4) imper- 
fect centering of the aperture and (5) PSF variations over 
the frame. These were minimised by using a sufficiently 
large aperture. To maximise the signal-to-noise ratio in 
aperture photometry, one usually chooses an aperture of 
size approximately equal to the H WHM (half-width at 
half- maximum) of the PSF (Howell 1989| ). It is not neces- 
sary to get all of the light - or apply an aperture correction 



- with relative photometry, provided that in a given frame 
a single-sized aperture is used and all stars have the same 
PSF. When there are PSF variations (as was the case for 
our data), a larger aperture must be used. This reduces 
the contribution from the informal errors (4 & 5) at the 
expense of increasing the formal errors (1~3). Thus the 
threshold for variability detections in any science star will 
be at a higher magnitude limit - but more reliable - than 
when using a smaller aperture. We found that an aperture 
size of about 4 x HWHM reduced the informal errors to 
well below the 0.5% level. Yet other sources of error (such 
as imperfect charge transfer efficiency and nonlinearity of 
the CCD flux sensitivity) are also below this level. 

After establishing reliable reference stars, the x^ test 
was applied to the science star. If it failed this test (gave 
evidence for variability), the next stage was to determine 
whether there was evidence for this variability being pe- 
riodic. For this purpose we evaluated the modified peri- 



odogram described by Scargle (1982). In contrast to the 



classical periodogram, it can be shown that the modifled 
periodogram is equivalent to a least-squares flt (in the 
time domain) of sinusoids to the data (Lomb 19761 ) . Fur- 
thermore, the method allows a simple determination of 
the significance of any peak in the periodogram based on 
the known noise in the data.[| We consider a peak to be 
significant only at the 99% level (i.e. when the probability, 
p, that a peak is due to noise is < 0.01). If a significant 
peak is found, the time series is phased to the detected 
period and examined to see if it shows cyclic variability. 



4. Results 

We now have four criteria which the data must satisfy 
before a science star is flagged as showing evidence for 
periodic variations: 

1. the science star has a significant Xs value 
(P(X^) < 0.01) 

2. the time series looks plausible 

3. the science star shows a significant (p < 0.01) peak in 
the periodogram 

4. the phased light curve at the detected period looks 
reasonable. 

These tests were applied to the six science stars in Ta- 
ble ^. All except 2M1145 failed at least one. From the 
results of the x^ test we can set approximate limits to 
the amount and timescale of intrinsic variability in these 
five stars at the epoch of the observations (Table l2|) . For 
each star in this table (except 2M1145) there is no signif- 
icant evidence (i.e. P{xi) is >0.01) for variability above 



Note tha t the Lomb-Scargle periodogram described by 



Press et al. (1992) differs from the periodogram of Scargle 



( 1982 ) by a constant factor. Thus the two methods will gen- 
erally give different significances for any peak in the peri- 
odogram. 
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Table 2. Approximate upper limits on intrinsic I band 
variability, /um is defined as the scatter of the RMS (root 
mean squared) values in the relative magnitude time series 
that would have been required before we obtained signif- 
icant (P(Xs) < 0.01) evidence for variability according to 
the x^ test. For each star the search was sensitive to vari- 
ations in approximately the range shown. 



Star 


Ilim 


period range 


2M0913 


0.070 


30 mins to 126.2 hours 


2M1145 


0.035 


30 mins to 74.8 hours 


2M1146 


0.025 


30 mins to 74.8 hours 


Roque 11 


0.045 


25 mins to 26.4 hours 


Teide 1 


0.050 


25 mins to 26.4 hours 


Calar 3 


0.050 


25 mins to 27.5 hours 



Ilim on the timescale shown. (This is not the same as say- 
ing that we are 99% confident that there is no variability 
above /um- Moreover, larger variability on timescales to 
which we were not sensitive could be present.) Strictly, it 
is only possible to put a single magnitude limit on a given 
range of periods if the uncertainties (e^.t) are the same 
for all t. This is not the case for our data due to changes 
in atmospheric conditions. Rather than defining different 
amplitude limits for different periods, we simply give an 
approximate limit for all periods. The upper limit period 
shown is the longest time base in our data. The short- 
est period is a sort of Nyquist period for the time series, 
and is slightly more than twice the shortest separation be- 
tween frames. We emphasise that these figures should not 
be viewed as strict magnitude/period sensitivity limits. 

The star 2M1145 passed all four tests for variability, 
and showed a significant {p < 10^^) peak in the peri- 
odogram at about 7 hours. The phased light curve is shown 
in Fig |l|. To be sure that the detected period was not due 
to a reference star, we calculated the periodogram of the 
relative magnitudes of each reference star (relative to the 
other reference stars, as described above). No reference 
star showed any even marginally significant peak. Addi- 
tional checks on the variability detection in 2M1145 were 
made by using only a subset of the good reference stars 
in the reference set, and using slightly different aperture 
sizes in the original photometry. In all cases significant 
variability was detected (according to the above four cri- 
teria) and the determined periods were the same to within 
1%. 

The RMS (root mean squared) scatter of the relative 
magnitudes in Fig |^ is 0.038 mags, and the amplitude of 
a least-squares fit sinusoid {A x sin[wi-|-0]) is 0.040 mags. 
This latter value may be a slight overestimate of the ampli- 
tude of intrinsic variability on account of noise: the larger 
the aperture, the more noise in each measurement, and so 
the more likely it is that a larger amplitude is observed. 
This does not mean, however, that the detection is just 
due to noise, as with a much larger aperture the time se- 
ries does not meet the variability criteria described above. 
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Fig. 1. Light curve of the (zero mean) relative magnitude 
time series for the L dwarf 2M1145 phased at the detected 
period of 7.12 hours. Two (identical) periods are shown. 
The RMS scatter of the variations is 0.038 mags, and the 
amplitude of the overplotted least-squares fit sinusoid is 
0.040 mags. 



This problem with least-squares fitting could be overcome 
using more robust techniques, but we choose to acquire 
a more extensive data set before making an improved de- 
termination. Nonetheless, Figl^ shows evidence of periodic 
variation beyond the size of the error bars. 



5. Discussion 

The most plausible explanation for the observed periodic 
variation in 2M1145 is a rotational modulation of the 
emitted flux. Assuming a radius of O.IM0 (Burrows et 



al. [1997D , and rigid rotation, the period of 7.12 hours im- 
plies an equatorial rotation velocity of 17kms~^. This falls 
in the range of rotation speeds for 93 field M dwarfs mea- 
sured by Delfosse et al. ( 1998| ) (all except one with wsini 
in the range < 2-32 kms~^), but is smaller than the range 
for 9 Pleiades M5-M6.5 dwarfs observed by Oppenheimer 



et al. (1997) (37 < wsini <65kms~^). However, our data 
does not give unambiguous evidence for rotation at this 
speed. All we can say for certain is that we have evidence 
for periodic variability which is not present in the refer- 
ence stars, and therefore is probably intrinsic to 2M1145. 
But given that we only have 28 points in our time se- 
ries spread over several periods, confirmation of this pe- 
riod is required with a more extensive data set, preferably 
with smaller error bars and in more than one filter to 
provide somewhat independent measurements of the pe- 
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riod. Observations should also be carried out over at least 
two (and preferably three or four) complete periods. We 
also stress that the period determination method assumes 
that the time series is stationary, in particular that the 
period and amplitude of the variations are constant: any 
evolution of surface features over the timescale of the ob- 
servations would interfere with the interpretation of the 
periodogram. Additionally, multiple surface features may 
not give rise to a single sinusoidal modulation (and indeed, 
other peaks were present in the periodogram). 

Given that the observations have only been carried out 
in a single filter, we can only speculate about the cause of 
the modulation in 2M1145. If its Ha emission (Kirkpatrick 
et al. 1999) can be taken as evidence of magnetic activity, 
then the modulation could be the result of magnetically 
induced star spots. However, as we have only observed 
three L dwarfs we cannot draw any conclusions about the 
correlation between Ha emission and rotation speed at the 
bottom of the main sequence, particularly as the ampli- 
tude limit on one of the targets (2M0913) is rather high. 
The observed modulations in 2M1145 could alternatively 
be the result of inhomogenous dust clouds rotating across 
the stellar disk. To distinguish between these two possibil- 
ities it will be necessary to re-observe in multiple filters (or 
with time resolved spectroscopy) to measure the change 
in Tes, and in filters sensitive to high dust opacity. 
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